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ABSTRACT 

The final stages of mass assembly of present-day massive galaxies are expected to occur through the accretion 
of multiple satellites. Cosmological simulations thus predict a high frequency of stellar streams resulting from 
this mass accretion around the massive galaxies in the Local Volume. Such tidal streams are difficult to observe, 
especially in dense cluster environments, where they are readily destroyed. We present an investigation into 
the origins of a series of interlaced narrow filamentary stellar structures, loops and plumes in the vicinity of the 
Virgo Cluster, edge-on spiral galaxy, NGC 4216 that were previously identified by the Blackbird Telescope. 
Using the deeper, higher-resolution and precisely calibrated optical CFHT/MegaCam images obtained as part of 
the Next Generation Virgo Cluster Survey (NGVS), we confirm the previously identified features and identify a 
few additional structures. The NGVS data allowed us to make a physical study of these low-surface brightness 
features and investigate their origin. The likely progenitors of the structures were identified as either already 
catalogued VCC dwarfs or newly discovered satellites caught in the act of being destroyed. They have the 
same g-i color index and likely contain similar stellar populations. The alignment of three dwarfs along an 
apparently single stream is intriguing, and we cannot totally exclude that these are second-generation dwarf 
galaxies being born inside the filament from the debris of an original dwarf. The observed complex structures, 
including in particular a stream apparently emanating from a satellite of a satellite, point to a high rate of 
ongoing dwarf destruction/accretion in the region of the Virgo Cluster where NGC 4216 is located. We discuss 
the age of the interactions and whether they occurred in a group that is just falling into the cluster and shows 
signs of so-called "pre-processing" before it gets affected by the cluster environment, or in a group which 
already ventured towards the central regions of Virgo Cluster In any case, compared to the other spiral galaxies 
in the Virgo Cluster, but also to those located in lower density environments, NGC 4216 seems to suffer an 
unusually high heavy bombardment. Further studies will be needed to determine whether, given the surface 
brightness limit of our survey, about 29 mag arcsec"^, the number of observed streams around that galaxy is 
as predicted by cosmological simulations or conversely, whether the possible lack of similar structures in other 
galaxies poses a challenge to the merger-based model of galaxy mass assembly. 

Subject headings: galaxies: clusters: individual (Virgo) - galaxies: dwarf - galaxies: evolution - galaxies: 
interactions 



L INTRODUCTION 

Dwarf galaxies play a key role i n ACDM-large scal e 
structure formation scenarios, e.g. iDekel & SilBdlgsel) : 
iNavarro et al.l (119961) . Indeed, according to this model, the 
dominant process of galactic growth is the accretion of small 

' Based on observations obtained with MegaPiime/MegaCam, a joint 
project of CFHT and CEA/DAPNIA, at the Canada-France-Hawaii Tele- 
scope (CFHT) which is operated by the National Research Council (NRC) 
of Canada, the Institut National des Sciences de I'Univers of the Centre Na- 
tional de la Recherche Scientifique of France and the University of Hawaii. 

t Email: sanjaya.paudel@cea.fr 



galaxies, and the large stellar halos found around massive 
galaxies such as our own Milky Way have been assembled 
thanks to the infall (and merging wi th) of many small dwarf 
galaxies dBuUock & Johnston! 120051) . Prior to their capture, 
dwarf galaxies are often disrupted through large tidal forces 
generated by their host galaxies, which form stellar stre ams 
and shells around their hosts (See review by lDucll20T2l and 
references therein). Although such fine structures have a life- 
time of a few Gyr in isolated environments, they are much 
shorter lived in dense cluster environments, where they in- 
teract with the cluster potential and are quickly dispersed 
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out ( lMihosll200l iTal et al] 120091: lAdains et alj|2012h . Dis- 
rupted satellites may also be contributing stars to the intra- 
cluster light, where the debris generated by numerous inter- 
actions over the lifetime of the cluster is obser ved as a vast 
diffuse halo around the most massive galaxies (iRudick et al.l 
[2009; Mihos_e_tj.l, 2005,; J anowiecki et al. 2010l) . On the other 
hand, the cores of tidally stripped galaxies of intermediate 
mass may become compact ellipticals (cEs), and the nuclei 
of tidally stripped dwarfs may become ultra-compact dwarf 
galax ies (ISasaki et all 120071: iForbes et al]l2003t iHuxor et al.l 
l20Tlh. The scena rio is commonly known as tidal threshing 
(iBekki et alj|2003l) . 

Because of their being stretched out, disrupted dwarf galax- 
ies experience a strong dimming of their surface bright- 
ness and become very difficult to detect. Until recently, 
very few candidate "disrupted dwarfs" were known outside 
of the Local Group. With the advent of sophisticated in- 
struments and the development of data analysis techniques 
optimized to detect low surface brightness ( LSB) struc- 
tures, the number of detections is grow i ng (iForbes et al.l 
200l l^totfaez-Delgado et al.1 120091 ifoTot iDuc et al.ll20Tlt 
Miskolczi et alj|201 U iKoch et alj|2012l) . However, most of 
the observed features - filaments, warps, plumes and shells - 
are left-overs from rather old collisions, as their progenitors 
are no longer visible or identifiable. The deep imaging sur- 
vey of nearby spirals carried out by iMartmez-Delgado et al.l 
( 120101) with robotic facilities such as the 0.5m Blackbird tele- 
scope show a number of cases of dwarf galaxies currently 
undergoing tidal disruption. These include a pair of dwarf 
satellites of the Virgo spiral galaxy NGC 4216 which are as- 
sociated with extended tidal tails, whose brightest parts, cat- 
aloged as independent dwarf galaxies in Virgo Cluster Cat- 
aloge (VCC: iBinggeli et al. 1985), are visible on shallower 
SDSS images iMiskolczi et al.ll20l il) . 

In this paper we revisit this system, using images obtained 
as part of the Next G eneration Virgo Cluster Survey (NGVS, 
iFerrarese et al.ll20T2l) . The excellent image quality of the sur- 
vey allowed us to disclose new objects, especially a couple 
of filaments and disrupted dwarfs, and to perform a surface 
photometry and color analysis. 

NGVS observations and data reduction are presented in 
Sect. ID The results, in particular the photometric properties 
of the newly discovered features, are detailed in Sect. |3] In 
Sect.m we investigate the origin of the streams, their age and 
fate. From a global view provided by the NGVS, we consider 
the overall uniqueness of the NGC 4216 system within the 
Virgo Cluster 

2. OBSERVATIONS AND DATA REDUCTION 

The observations were carried out as part of the NGVS 
deep imaging survey using the MegaCam instrument 
(iBoulade et al. 2003) on the 3.6m Canada-France-Hawaii 
Telescope (CFHT). The observational strategy, data reduc- 
tion procedure and data quality control are described in 
IFerrarese et al.l (12012b . In brief, the multi-band NGVS survey 
benefits from the very wide field of view of the camera (~ 1° 
X 1°), and its very good image quality (0.6"-0.9") obtained 
through optimizing the observing conditions. The images are 
processed with a dedicated data reduction pipeline, Elixir- 
LSB, which is explicitly designed to optimize the detection 
of very low surface brightness structures such as intra-cluster 
ligh t (ICL), faint stellar streams and ultra-faint dwarf galax- 
ies dMagnier & Cuillandrel 120041 Cuillandre et al., in prepa- 
ration). The output products of Elixir-LSB are single frames 



TABLE 1 

Basic properties of NGC 4216 and its companion galaxies. 
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Note. — 


The coordinates, heliocentric radial velocities 
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phologies are from NED. The absolute B-band magnitudes are derived from 
the B-band magnitudes listed in the Goldmine databasfl An average dis- 
tance for the Virgo Cluster of 16.5 Mpc (corresponding to 80 pc arcsec"') 
was assumed, and d is the projected distance on the sicy from the center of 
NGC 4216. 
"http://goldmine.mib.infn.it 

and stacks in which the s ky background has been modeled 
and removed (lGwynll2008l: Gwyn et al. in preparation). Fur- 
ther image analysis has been done using standard routines of 
IRAlQ and IDL. 

For the study, we made use of NGVS images obtained in 
the g', i' and z' bandfl Total exposure times were 0.88, 0.57 
and 1 .2 hours respectively. Since the ^-band image is the most 
sensitive and the cleanest, i.e., it contains the least number of 
artifacts such as CCD imprints and reflection halos (see be- 
low), it was preferentially used for our photometric analysis. 

Surface photometry of very low surface brightness objects 
is notoriously difficult. The precision of their measured pa- 
rameters is affected by several factors, including the accuracy 
of the photometric calibration, sky background variation and 
artifacts due to imperfect optics. The photometric calibra- 
tion of the NGVS has been done using the SDSS photome- 
try with a level of accuracy corres ponding to a 1 % vari ation 
in the photometric zeropoint (see IFerrarese et all 120 121) . At 
the surface brightness limit of the survey (29 mag arcsec in 
the g-band), the main limitations to the LSB structure analy- 
sis and photometry are due to the reflection halos around the 
bright sources (stars but also the nuclei of bright galaxies with 
steeply rising surface brightness profiles), which show up as 
large disks that cannot easily be removed from the images 
(see an example of one of these halos in the upper-left corner 
of Figure [Til. We therefore exclude in our analysis the con- 
taminated areas where reflected light becomes dominant over 
the actual LSB feature. 

Foreground stars and compact background galaxies were 
removed applying a ring filter to the images, and residuals 
were manually subtracted with the IRAF task imedit. The 
apertures used to carry out the photometry of the faint fila- 
ments and plumes, shown in Figure|2] were defined visually as 
their surface brightness is too low for an automatic detection. 
The errors in the photometric measurements are to a large ex- 
tent determined by those in the sky determination. Although a 
master sky background is subtracted from each image, for our 
aperture photometry measurements we sampled the sky back- 
ground at multiple positions around the targets. The number 
of sky probes ranged between 5 and 100, depending on the 
size of the objects. Each individual sky region consisted of 
20 X 20 pixel boxes. The median value and errors were then 

^ Image Reduction & Analysis Facility Software distributed by National 
Optical Astronomy Observatories, which are operated by the Association of 
Universities for Research in Astronomy, Inc., under co-operative agreement 
with the National Science Foundation 

* noted as g, i, z in the rest of the paper 
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computed. 

3. RESULTS 

3.1. NGC 4216 and its environment 

Figure [T] shows a composite Meg aCam image of the field 
around NGC 4216 (cf. Fig. Ic in iMartmez-Delgado et al.l 
dMl)- The regions with the highest surface brightness 
are displayed in "true color" (combination of g, i and z 
band images) while for the low surface brightness regions, a 
monochromatic g-band image is shown. 

NGC 4216 is a barred spir al galaxy - it is classified as 
SAB(s)b in the RC3 cataloge dde Vaucouleurs et al.l [19911) - 
located in the outskirts of the Virgo Cluster at an angular dis- 
tance of ~4 degrees from the cluster center, MSTQ. It has 
an edge-on orientation with an inclination angle / = 85°. As 
shown in Figure [T] it s central regi ons are characterized by a 
prominent dust lane ( IChung et al.l l2009) and a red and com- 
pact bright nucleus. The galaxy is surrounded by a very dif- 
fuse stellar halo, and complex filamentary structures, includ- 
ing streams and plumes. The investigation of their nature is 
the focus of this study. 

The Hi map of NGC 4216 from the VLA (IChung et al.1 
120091) indicates an Hi extent which is somewhat smaller than 
in the optical: the ratio of the Hi and optical B-band radii, 
measured respectively at the 1 M© pc"^ Hi surface density 
level and the 25 mag arcsec"^ isophotal level in the B-band, is 

^ = 0.76. 

Ob 

The low radial velocity of NGC 4216 (131 kms"', which is 
significantly offset relative to the ^1300 kms"' mean veloc- 
ity of the cluster itself) indicates that it is moving relatively 
fast with respect to the cluster center It is the most mas- 
sive galaxy within a circle of 400 kpc projected radius; the 
nearest galaxy of higher mass is M 99 at a projected distance 
of 420 kpc. Within that field of view, there are 2 early-type 
galaxies, and 5 late-type star-forming galaxies with radial ve- 
locity differences less than ±300 kms"' with respect to that of 
NGC 4216. The basic properties of NGC 4216 and its clos- 
est companions are listed in Table [T] In the close vicinity of 
NGC 4216, towards the North lies VCC 165, at a projected 
distance of 19.2 kpc. VCC 165 is classified bv lBinggeli et alj 
(II985') as an SOi galaxy. Another companion is located fur- 
ther (56.5 kpc) to the North East, at the same radial velocity as 
VCC 165: the spiral galaxy NGC 4222. On our deep images, 
the stellar halo of the galaxy NGC 4222 overlaps with the re- 
flection halo of a bright star. South-East of NGC 4216, at 
58.0 kpc distance from its centre, at a similar radial velocity, 
lies a nucleated early type dwarf galaxy VCC 200 (Cote et aLj 
l2006h . All these companions are identified on Figure |2] 

3.2. The streams around NGC 4216 

Figure |2] identifies the main streams and remnants of dwarf 
galaxies in the vicinity of NGC 4216 and VCC 165. The 
NGVS images clearl y show the main long filament F l, which 
was first detected bv lMartmez-Delgado et al.l (120101) . We also 
se e the fainter structures present on the deep optical image 
of IMartmez-Delgado et al.l ( 120101) such as the bend in Fl to- 
wards its Easter tip, the long narrow filament to the North of 
the system (F2) and the plume-like filaments that stick out 
East of NGC 4216 (F5 and F6). Besides these previously 

^ The virial radius of Virgo's A subcluster, centred on M87, is 5.4° 
(McLaughlin 1999). 



known structures, new features are detected, such as a fila- 
ment apparently wrapping around VCC 165 (F4), a loop East 
of NGC 4216 (F3) and a possible bridge between NGC 4216 
and VCC 200 (F7). 

With respect to previous surveys, the gain in spatial res- 
olution provided by CFHT/ MegaCam allows us to identify 
the dwarf galaxies presumably at the origin of these fila- 
ments. Finally, the photometric calibration precision of the 
multi-band NGVS allowed us to determine the photometric 
properties and colors of these various features, which so far 
had not been possible. The photometric data on the streams 
and progenitors are presented respectively in Table 0& |3] 

The most prominent tidal stream, Fl, East of NGC 4216, 
has an angular extent of ~12.5 arcmin on the sky, which cor- 
responds to ~ 60 kpc at the assu med mean distance of the clus- 
ter, 16.5 Mpc (iMei et al.ll2007h . The filament starts from the 
South-East corner of NGC 4216 and near its easternmost tip 
it makes a sudden turn towards the West-North with an an- 
gle of about 120 degree. Although further away it gets lost 
within the halo of a bright red star, when it emerges again, 
it seems to join the filament F2, the second longest filament 
in the field. North of NGC 4216. Unfortunately, the stellai- 
halo glare in the region between the two filaments does not 
allow us to say with certainty that a physical connection ex- 
ists. No such bridge is visible on slightly shallower image of 
IMartmez-Delgado et al.l (120101) , where the stellar halo is less 
obtrusive. The width of the stellar stream Fl appears constant 
along most of its length, about ~3 kpc. It broadens after the 
bend. As shown in Figure |4] the average surface brightness 
of the filament (excluding regions close to the dwarf galaxies 
along it) is slightly below 28 mag arcsec"-. 

The NGVS image discloses another rather broad stream, 
F3, with a surface brightness below 28.5 mag arcsec"^ (See 
Figure |3]l, which sticks out of NGC 4216, then bends to cross 
filament Fl and finally falls back towards the host galaxy, near 
F7. 

The very broad filament F7 points in the direction of the 
dwarf companion VCC 200, but a physical connection be- 
tween the two galaxies cannot be firmly established. The 
plume-like stellar structures which apparently stick out from 
the main body of NGC 4216 on its western side (F5 and F6) 
may in fact be filaments which connect to F1-F3 behind the 
galaxy. They have a nearly uniform surface brightness of 27.5 
mag arcsec"-. 

Other filamentary structures can be seen (best on Figure ^ 
around VCC 165, the closest companion of NGC 4216. On 
the displayed g-band image, a model of VCC 165 obtained 
with the IRAF task bmodel has been subtracted to further 
enhance the filaments. The most prominent one, the 6-kpc 
long structure labeled F4, wraps around the dwarf companion 
on one side, and points towards the galaxy center on the 
other There is a further hint that the galaxy extends further 
towards the East in the direction of the spiral. VCC 165 
itself, despite its apparent proximity to NGC 4216, does not 
exhibit strong signs of morphological perturbation. The map 
obtained subtracting the ellipse model (see Figure|6) does not 
show prominent residuals in its outer regions. In the inner 
regions, spiral-like features (or narrow dust lanes) may be 
observed. 

We measured the g - i colors for the most prominent fea- 
tures. The /-band photometric measurements were particu- 
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Fig. 1 . — Composite NGVS image of the field around NGC 4216. A monochromatic g-band image for which the faintest point-like stars have been subtracted 
is shown with a grey scale. For regions above a surface brightness level ~24 mag arcsec"^, and in empty sky regions far from the main galaxy, a true color image 
(composite of g, i, z-band images) is superimposed. North is up. East is left and the field of view is 20 X 17 arcmin. The image in the inset on the upper right 
comer is a zoom towards VCC 165, the closest companion of NGC 4216. 



larly delicate as this band shows several residuals at low sur- 
face brightness levels, such as chip gaps. The contamination 
by stellar halos is also more severe in this band. Even the nu- 
cleus of NGC 4216 generates an artificial, disk-like structure 
with a radius of nearly 3 arcmin, which is almost invisible in 
the g-band, but becomes noticeable in the / and z bands. Re- 
gions affected by the stellar or galactic reflection halos were 
excluded from our analysis. The colors are listed in Table |2] 
The g-i color for Fl is uniform along its length and has an 
average value g-i = 0.9 ±0.1 mag. Interestingly, this value 
agrees, within the errors, with that of filament F2, possibly 
suggesting a common physical origin for the two filaments. 
On the other hand, there is a significant difference of 0.4 mag 
between filaments F4 and F2, although both structures touch 
each other on the projected image. As argued later, they may 
have different progenitors, and indeed their different orienta- 



tion supports this hypothesis. 

3.3. Stellar cores within the filaments 

Several objects of higher surface brightness are visible to- 
wards the streams: such cores are likely the remnants of 
tidally disrupted dwarf galaxies and as such appear as good 
progenitor candidates of the filaments' stellar material. 

As shown in Figures |2] and |3] the most extended filament, 
Fl, hosts three relatively compact stellar objects. The most 
luminous object. A, is also the most extended (see Figure O. 
It has a ~24.5 mag arcsec"- central ^-band surface bright- 
ness, which in principle is high enough to have it included in 
the Virgo Cluster Catalog (VCC 197). Modeling the galaxy 
with ellipses (using the IRAF task ellipse and bmodel), we 
found that it is stretched out along filament Fl and contains 
a bar-like structure along the minor axis, and a compact nu- 
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TABLE 2 

Properties of the filaments (F1-F7) and external halo regions 
(G1-G2) OF NGC 4216. 



Fig. 2. — Identification of tlie principal companion galaxies, stellai' streams, plumes and progenitor candidates around NGC 4216. The apertures used to cany 
out the photometry are delineated on the displayed g-band image. Note that for objects A, B and C, the boundaries correspond to a surface brightness of 27.3 
mag arcsec"^ , the average surface brightness of filament F 1 . G 1 and G2 are two selected regions in the extended stellar halo of NGC 42 1 6. The dashed line near 
the extended stellar halo represents a possible, but not yet proven, connection between Fl and F2. 

cleus. Modeling t he galaxy with a Sersic model with GALFIT 
(iPeng et alj|200l . we obtained a best fit with an exponential 
profile and a very large effective radius of ~5 kpc, a conse- 
quence of the tidal threshing. 

Object B - also in the VCC (VCC 185) - with an effective 
radius of 0.4 kpc and an ellipticity of 0.29, is more compact 
and rounder Although located at the crossing of filaments Fl 
and F3, its major axis is aligned with neither of these streams. 
Finally, object C is an uncataloged faint dwarf located near 
the base of the tail. It has a round shape and does not show 
any clear evidence of a tidal interaction. The integrated g-i 
colors of both A and B are in very good agreement with that 
of filament Fl, while C is 0.2 mag bluer. Given their rather 
red colors, none of the galaxies show evidence of containing 
young stars. A is not detected on a deep 2600 sec archival 
GALEX/NUV image jBoselli etalj|2011b . while B shows 
very weak UV emission . At Arecibo, the ALFALFA Hi sur- 
vey jHavnes et al.ll201 lb did not detect any atomic hydrogen 
emission towards the three passive dwarfs. The much deeper 
AGES fields did not cover the NGC 4216 area (Tavloret a!] 
l20T2allbh . 

The brightest condensation along Filament F4, West of 
VCC 165, turns out to be a chance superposition of a back- 
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Note. — The color g — i was determined from aperture photometric mea- 
surements in areas not contaminated by instrumental artifacts. < > is the 
average surface brightness within the areas defined in Figure|2] excluding the 
progenitors (i.e. A, B & C). The associated en'ors are computed from the stan- 
dard deviation of the individual sky background measures. The absolute blue 
magnitude in the B band, Mb, was derived from the extrapolated observed 
magnitude, i.e., the average surface brightness times the area, converting the 
MegaCam magnitudes to the standard Johnson system using the calibration 
of Joidi et al. (2006) and assuming a distance of 16.5 Mpc iMei et al. 200^. 
The stellar masses, M«, were estimated fro m the extrapolate d absolute mag- 
nitudes assuming the mass to light ratios of lBell et alj l l2003h . 
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Fig. 3. — g-band surface brightness map of the field around NGC 4216. The surface brightness scale in mag arcsec"^ is shown to the right and the field of view 
is the same as for Figure |2] The white contour corresponds to the faintest level of 21 -cm Hi line emission in NGC 4216 detected with the VLA bv lChung et alj 
<20Q9l) . 




30 

Length, kpc 

Fig. 4. — Variation of the median g-band suiface brightness along the length of filament Fl, from dwarf galaxy C near' NGC 4216 to its easternmost tip where it 
gets lost within the halo of star. Photometric measurements were obtained collapsing the filament along its width and median-combining the data. The positions 
of the dwarf galaxies A, B and C (see Fig. 2) are indicated. 



ground low-mass AGN at a redshift of 0.043, according to its 
SDSS spectrumQ. The second brightest object in that region 
is an extended structure that lies in the continuation of F4. 



^ The AGN, SDSS J121551.26-H131303.4, belongs to a rich large-scale 
filament of galaxies. Its emission lines classify it as a Seyfert 2. With an 
absolute g-band magnitude of -18.2, it is among the least massive AGNs 
found in the whole SDSS survey. A visual inspection of all spectroscopically 
confirmed SDSS AGNs below z < 0.05, less luminous than Mg > -19 mag, 
did not reveal other objects with tidal tails apparently associated to them. 



Given its much higher surface brightness than that of filament 
(see Figures|5]and|7]i, we argue here that it could be in fact the 
main-body remnant of an almost completely disrupted dwarf, 
which we have labeled on the figures as D. The fact that D 
and F4 have similar colors is consistent with this hypothesis. 
A small compact central core, labeled as N, can be seen to- 
wards D; it has the same color and could be the nucleus of 
the latter The B-band absolute magnitude of D - within the 
boundaries defined in Figure|6]- is -11.3 mag, i.e. in between 
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Fig. 5. — g-band images of the dismpted dwarfs A, B, C and D. For objects 
A, B and C, located along filament Fl, the images have been rotated such 
that the local orientation of the filament is parallel to the horizontal axis. 
D is a disrupted dwarf' galaxy candidate located West of VCC165. On the 
displayed image, the latter galaxy has been subtracted using an ellipse model. 
All images are displayed with the same intensity and spatial scales. 

that of dwarfs A and B. 

The mass of the stars still bound to the dwarfs (A+B+C+D) 
amounts to less than one tenth of the total stellar mass of the 
filaments (Fl-7). Assuming that all the stellar material in fila- 
ment Fl was pulled out from galaxy A, then the galaxy must 
originally have had an absolute magnitude of = -15.6 
mag, a value in the r ange of bright earl y-type dwarf galaxies 
in the Virgo Cluster (iLisker et al.ll2007h . 

4. DISCUSSION 

In the Carnegie Atlas of Galaxies (ISandage & Bedkell994h . 
the spiral galaxy NGC 4216 is described as "one of the most 
famous galaxies in the sky because it is often used in text- 
books to illustrate the bulge, disk, dust content, and spiral 
pattern of typical luminous galaxies near the middle of the 
Sb classification sequence". The deep NGVS images we have 
presented here show many additional features around this pro- 
totypical galaxy, in particular a complex network of interlaced 
plumes, streams and filaments. Although a detailed interpre- 
tation is not straightforward, the large number of low sur- 
face brightness stellar tails is consistent with predictions from 
numerical simulation s about the mass assembly of galaxies 
jJohnston et al.ll2001l: iBuTlock & Johnstonll2005l) . This raises 
the question why all massive galaxies do not show a similar 



TABLE 3 

Properties of the dwarf galaxies identified along the filaments. 



Object 


Mb 




e 


log(M*) 


Re 




mag 


mag 




Mq 


kpc 


A(VCC 197) 


-12.1 


0.88 ± 0.02 


0.79 


7.2 


4.5 


B(VCC 185) 


-10.2 


0.86 ± 0.02 


0.29 


6.7 


0.4 


C 


-09.6 


0.62 ± 0.04 




6.0 


0.3 


D 


-11.3 


1.05 ±0.02 




7.1 




(N) 


-07.0 


1.05 ±0.04 


0.63 


5.3 




VCC 165 


-15.5 


1.22 ±0.01 


0.15 


8.9 


1.3 



Note. — The absolute blue magnitude and colors were detennined from 
photometric measurements made within the apertures defined in Figure |2] 
The ellipticities e were derived using the IRAF ellipse task. The effective 
radius are derived with a 2D fitting of the light profile with a Sersic model 
made with GALFIT (Peng et al. 2002). D is the disrupted dwarf galaxy near 
VCC 165 (see Fig.|6), while N is its candidate nucleus. 




Fig. 6. — g-band image of the field around VCC 165. The main features 
have been labeled: the streams F2 and F4, the progenitor candidate of F4, the 
disrupted dwarf D. with its possible nucleus, N, and the background AGN. 
See the inset in Figure[T]for another view in color on this complex structure. 




2 4 6 8 

Length, kpc 

Fig. 7. — Surface brightness profile along filament F4. The black solid 
curved line in Fig|6]dehneates the region traced for the profile. 
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level of fine-structures. In the following, we propose scenar- 
ios accounting for the various features seen around NGC 42 1 6 
and discuss what is so special (or not) about this galaxy. 

4.1. Real structures or artifacts 

Disentangling real features from artifacts is essential to the 
goals of this paper. Although multiple internal reflections in 
CCDs and stellar halos also produce low surface brightness 
structures, these instrumental features usually have somewhat 
regular geometric shapes, like disks, long thin straight lines 
or grids, while the shapes of the cosmic structures of inter- 
est to us are significantly different in appearance: broad fila- 
ments, arcs and irregular plumes. In the case of NGC 4216, a 
number of the faint struct ures we note d are (barely) visible on 
SDSS images (Miskolcz i et alj|201 II) and more convincingly 
detected bv iMartmez-Delgado et al.l (120101) . 

As a reminder, the principal features observed around 
NGC 4216 and possibly associated with collisional debris are 
the following: 

(a) two long filaments (Fl and F2), with frail evidence that 
they are physically connected, forming a possible loop in the 
North- West quadrant of NGC 4216. Three dwarf galaxies (A, 
B and C) are found projected onto, and we argue, associated 
with filament Fl; 

(b) a less extended very low surface brightness arc-like struc- 
ture east of NGC4216 (F3). The dwarf sateflite B is found at 
the intersection of F3 and Fl; 

(c) a relatively high surface brightness filament wrapping 
around the dwarf galaxy companion VCC 165, likely ema- 
nating from a newly discovered disrupted dwarf (D); 

(d) several broader and short filaments coming out from main 
body of NGC 4216 (F5-7). 

4.2. Origin of the filamentary structures 

Undoubtedly, all the stellar filamentary structures observed 
in the field of NGC 4216 have been shaped by tidal forces. 
However the origin of their stars may be debated. Given the 
high frequency of fine structures around the spiral, one may 
wonder if they were all produced by a single event: either an 
old merger with a relatively high mass companion, or ongo- 
ing tidal interaction with dwarf galaxies. In these two cases, 
the dwarfs seen along filament Fl would not be the original 
progenitors of its stars, but they would instead be tidal dwarf 
galaxies (TDGs) bom in situ from the collisional debris. In 
fact, the three aligned dwarfs A, B and C show similarities 
with the three also aligned TDG candidates recently discov- 
ered along a prominent tidal tail associated with the elliptical 
galaxy NGC 55 57, which was interpreted as an old merger 
jDuc et ani20I Ih . However the latter are gas-rich, contrary to 
the NGC 4216 dwarfs. 

NGC 4216 has at least one relatively massive companion, 
the edge-on spiral NGC 4222, and two luminous dwarf satel- 
lites, VCC 165 and VCC 200 (see Table The latter two 
galaxies are likely not massive enough to tidally disrupt their 
host (though the filament F7 pointing towards VCC 200 is 
intriguing). An on-going collision between the two spirals 
would generate tidal tails in t he pl ane of their stellar disks 
(see review by [Due & Renaud ,2013h . In fact, all the new fea- 
tures discovered around the galaxy (the stellar protuberances, 
plumes and filaments) are located along its minor axis, i.e., 
perpendicular to its main plane, and not towards the outermost 
regions of its stellar disk plane, where the isophotes remain 
regular (see in particular Figure [T). In Hi. .Chung et al..(.2009i) 



did not find any clear indications of interactions between the 
two galaxies either, though they noted mild distortions in the 
kinematics at the edge of the disk in both systems. A ma- 
jor merger in NGC 4216 is thus very unlikely. Furthermore, 
the outer regions along the major axis of the galaxy, Gl and 
G2 (see Figure |2]i, are slightly bluer in ^ - / by 0.2 mag than 
the regions along the intergalactic stellar streams. This likely 
means that the stellar streams and outer disk of the galaxy do 
not share the same stellar populations, and thus that the stars 
in the filaments did not originate from the spiral. 

Alternatively, the filamentary structures may result from 
dwarf galaxies tidally stirred by NGC 4216. In that case, the 
tails are shaped by both tid al forces and the orbit of the satel- 
lite around the host galaxy ( IVarghese et al.ll20TTl) . Inside sev- 
eral of the filaments, stellar bodies with higher surface bright- 
ness are found that are good filament progenitor candidates. 
The most promising one is object A, which shares with fila- 
ment Fl its orientation and color. Scrutinizing Figure [T] one 
may see that the filament has a S-shaped structure close to 
the dwarf, due to the presence of two slightly misaligned tails 
on each side of the galaxy, possibly the leading and trailing 
tidal tails. This provides an additional evidence that A is the 
progenitor. 

On the other hand, it is hard to determine whether object 
B belongs to filament Fl or F3, or if it is just a chance 
superimposition. Its main body is not aligned with any of the 
filaments while its g-i color agrees with that of filament Fl 
(the S/N of F3 is too low to determine its color). If objects 
B, as well as C, are indeed associated with Fl, their presence 
along a single filament remains mysterious. It is unlikely that 
their common progenitor broke in two or even three stellar 
cores. Alternatively, if A, B and C were all pre-existing 
objects, and subsequently tidally disrupted by NGC 4216, 
one would a priori expect them to be originally on different 
orbits and therefore to generate filaments with different 
orientations. Note however that according to cosmological 
models, satellite accretion occurs along privileged directions. 
The dwarfs might thus belong to a common cosmological 
stream. A chance superimposition of B and C is obviously an 
hypothesis worth exploring. Having spectroscopic data, and 
in particular information on the radial velocity and chemical 
abundances, would help to determine the nature of the dwarfs. 
Alternatively, B and C might be substructures born in situ 
in the tidal tail of A. Such a formation mechanism would 
differ from the one at the origin of TDGs, as it would not 
involve a major wet merger. It would rather resemble the one 
responsible for the formation of the stellar clumps that are 
commonly fou nd in the tidal streams of MiUcy-W ay Globular 
Clusters (e.g.. iMastrobuono-Battisti et al.ll20I2l) . Note that 
in our own local group, the distribution of dwarf satellites 
inside narrow di sks, t he so-called disks of satellites, (e.g. 
[Pawlowski et al. [20121 ; ITbata et al. 2013), is also puzzling 
and no satisfactory solution to account for them has yet been 
found. 

Finally, the origin of the filamentary structures that do not 
show any obvious progenitor, such as F2 (unless it is con- 
nected to Fl and shares the same progenitor), and of the three 
plumes-like filaments (F5-7) is even less constrained. The 
progenitors might have been totally disrupted or their rem- 
nants might already have been swallowed by the massive spi- 
ral. In fact, object D near VCC 165 is likely a case of a dwarf 
in the process of total disruption. Whether it is disrupted fol- 
lowing its interaction with VCC 165 (around which its tidal 



Destruction of dwarf galaxies around NGC 4216 



9 



tail, F4, wraps) or with the more massive but slightly more 
distant galaxy NGC 4216, or both, is unclear VCC 165 itself 
does not seem to be affected by gravitational interactions. 

4.3. An unusual on-going bombardment in a cluster of 
galaxies 

The discovery of numerous filamentary structures around 
NGC 4216, and tidally stretched stellar clumps within them, 
revealed an impressive number of in-falling satellites caught 
it the act of disruption. When has such a bombardment started 
and where in the cluster did it occur? 

4.3.1. A sub-structure within the Virgo Cluster? 

Could NGC 4216 belong to a sub-structure, like a group 
that has been accreted by the Virgo Cluster? Within the clus- 
ter, NGC 4216 is located in a moderately dense environment: 
the local projected number density of galaxies is nearly half 
that of the cluster core. The spiral is the most luminous galaxy 
within this sub-structure. For each Virgo galaxy more lumi- 
nous than Mb = -18, we estimated in the SDSS dataset the 
number of companions within a 100 kpc radius and radial ve- 
locities of ±300 kms~'. NGC 4216 turns out to have a signif- 
icantly higher number of faint companions (M^ > -18, with 
spectra available), compared to its fellow galaxies. The prob- 
ability of an increased number of dwarfs being demolished is 
thus higher there. The statistics on the ratio of dwarfs to giants 
is still to be extended to the many faint and ultra-faint dwarfs 
that are being discovered by the NGVS survey. Low veloc- 
ity collisions able to generate multiple tidal tails are expected 
in groups. As a matter of fact, examples of tidal interac- 
tion and possibly ram-pressure have been observed in several 
in-falling groups of cl usters, which are su ffering a so-called 
"pre-processing" (e.g.. lCortese et al1l2006l) . during which the 
properties of their galaxies are changed within the group be- 
fore being affected by the cluster environment. 

4.3.2. Time scales for the satellite accretion/destruction 

A time scale for the satellite accretion/destruction process 
may be obtained from the survival time of the debris: whereas 
around isolated systems, streams may remain visible for sev- 
eral Gyrs, within a cluster, the life expectancy of tidal de- 
bris is largely reduced due to dynamical heating and mix- 
ing (iRudick et alj 120091) . They should evaporate on a time 
scale shorter than the clu ster crossing time, typically 1 Gyr 
jTrentham & Tullvl 120021) . This suggests that the bombard- 
ment has either started within the Virgo Cluster itself or in 
a pre-processed group that was accreted by the cluster only 
recently. 

The gas content of the NGC 4216 system somehow dif- 
fers from that of the pre-processed groups so far reported 
in the literature. The filamentary structures detected in 
claimed pre-processed groups (e.g., ICortese et al.l 120061) are 
usually gas-rich, host Hll regions, and are thus probably star- 
forming. There is no evidence for intergalactic Hi gas around 
NGC 4216, at least in the ALFALFA single-dish surve>Q 

The tails are red, presumably made of old stars. The dwarf 
satellites around NGC 4216 had their gas stripped long ago, 

' Considering that some large-scale gas features of low col- 
umn density such as the Virgo Hi plume discovered by the WSRT 
jOosterloo & van Go rkom 2005) have been missed in single dish surveys, 
deeper imaging studies with interferometers which will not resolve out dif- 
fuse gas will help to ensure the absence of intergalactic gas. 



either due to long-lasting interactions with the massive galax- 
ies in the group, or because globally the group members had 
their gas reservoir stripped before the accretion event. As 
noted in Sect. the Hi disk of NGC 4216 seems trun- 
cated, which is normally interpreted as a consequence of 
ram pressure stripping in the cluster environment: the outer 
gaseous halo of a spiral galaxy is stripped away when it 
passes t hrough a high densi t y medium such as the co re of a 
cluster dVollmeret all 120011: iBoselU & Gavazzil 120061) . This 
might indicate that the NGC 4216 group already suffered a 
ram-pressure stripping episode, and therefore that is has al- 
ready been in the cluster for a long time. The ram pressure 
needed to strip a spiral down to the observed outer Hi ra- 
dius Rhi {15 kpc) is ^ J:hiVI,/R = PicmV^^,, where Eh/ is 
the Hi surface density (3 Mq/ pc^), ¥,■„, the rotation veloc- 
ity (250 km s"') and Vgai the galaxy velocity in Virgo cluster 
(1500 km s"'). Thi s yields an ICM density, picM of 3 x 10"^ 
cm"^. According to lSchindler et al.l (119991) . such a density oc- 
curs at a cluster radius of about 0.5 Mpc or 1.7 degrees. At 
the present epoch, NGC 4216 is located at projected distance 
of 2.5 and 3.6 degrees respectively from the central cluster 
galaxies, M86 and M87. Thus, the value of its Hi truncation 
radius suggests a stripping episode that occurred at least 1 Gyr 
ago, given the galaxy velocity, and closer to the cluster core, 
which would imply that the galaxy has ventured through that 
environment in the pasfl The bombardment of NGC 4216 
and the multiple dwarf destruction must then have occurred 
later on. 

Alternatively, the Hi truncation of NGC 4216 might be 
possible without invoking the Virgo Intra Cluster Medium, 
if the group were sufficiently massive and X-ray luminous 
dSengupta et al.ll2007l:lkern et alj|2008l) . There is however no 
evidence for extended X-ray emission at the location of the 
galaxy to support the hypothesis of an i n situ ram pressure 
event (See Figure 7 in IChung et al.ll2009h . Finally one should 
note that one of the main characteristics of the gas property of 
NGC 4216, more striking than the disk truncation, is its over- 
all low Hi surface density for a spiral galaxy of its size. Such 
a global depletion of Hi mi ght in fact be due to the effect of 
turbulent viscous stripping (INulsenll 19821: IChung et aDl2009l) 
rather than ram pressure stripping. 

Taking furthermore into account the fact that the companion 
galaxy, NGC 4222, has a regular Hi distribution, it is still 
difficult to firmly conclude that the group of NGC 4216 did 
experience a ram pressure event during a past journey through 
or near the Cluster Core. 

4.3.3. How unique is NGC 4216? 

How special is NGC 4216 within the Virgo Cluster? As 
part of the NGVS project, we have identified in a systematic 
way low-surface brightness stellar tails, and cases of tidally 
disrupted dwarfs. Our preliminary analysis over most of the 
Virgo Cluster area indicates that there are only a rather small 
number of them, especially around spiral galaxies (Due et al., 
2013 in prep.). This may be an additional argument in agree- 
ment with a recent accretion of NGC 4216 as part of a group. 

The true distance and location of NGC 4216 within the Virgo Clus- 
ter is not known. An estimate of the distance based on Surface Brightness 
Fluctuations is only a vailable for one possible group member, VCC 200. 
IBlakeslee et alj 120091) determined a distance of 18.2 ±0.6 Mpc for this 
galaxy. The group thus may be located further from M87, at ~2 Mpc (cal- 
culated 3D distance using radial and projected distance 1.7 and 1.2 Mpc re- 
spectively). This yields a travel time from the centre of Virgo Cluster to the 
current position of nearly 2 Gyr. 
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On the other hand, even in the field and in groups, a bom- 
bardment like the one suffered by NGC 4216 seems pretty 
rare. N one of the spirals mapped by iMartmez-Delgado et al.l 
(1201 Oh have a similar number of filaments and disrupted 
dwarfs around them. In some systems, what appears to be 
multiple loops may in fact be due the wrapping of a sin- 
gle tid al stream around its host (e.g.. lMartmez-Delgado et al] 
120091) . In the case of NGC 4216, several independent fil- 
aments are observed; they emanate from at least 2, likely 
4, different disrupted dwarfs. The on-going deep imaging 
surve y with MegaCam of Atlas3D field galaxies (Due et al] 
1201 II) . which has a similar depth and sensitivity to low-surface 
brightness structures as the NGVS, has so far also failed to 
detect the same wealth of debris of minor mergers observed 
around NGC 42 1 6 . On the other hand, multiple stellar streams 
and dwarf satellite proge nitors have been identified in the halo 
of our own Milky Way (llbata et al.ll200lt lYannv et al.ll200l 
iGrillmair 2006) . They are however of much lower surface 
brightness. Were the NGVS able to reach the surface bright- 
ness limit achievable with resolved stellar populations - 32 
instead of 29 mag arcsec"^ - then the number of filaments and 
satellites is expected to strongly increase according to numer- 
ical simulations (iBullock & Johnstoa2005i, Michel-Dansac et 
al., 2013, in prep.). 



5. CONCLUSIONS 

We have presented a study of the properties and origins of a 
complex series of interlaced narrow filamentary stellar struc- 
tures, loops and plumes located in the vicinity of the edge-on 
Virgo Cluster spiral galaxy NGC 4216. They were detected 
on multi-band extremely deep optical Next Virgo Cluster 
Survey (NGVS) images. Some streams had alread y been 
identified in previous surveys, e.g. Martmez-Delga do et al.l 
(l2OT0h . Several others, with g-band surface brightness levels 
fainter than 28.5 mag arcsec"- are reported here for the first 
time. 

Whereas the tidal origin of the filamentary structures seems 
rather obvious, determining the progenitors of their stars is 
less straightforward. The high resolution of the MegaCam 
camera, besides its high sensitivity to low surface brightness 
features, allowed us to identify progenitor candidates and to 
further explore the complexity of the system. This is illus- 
trated by the detection of satellites of satellites: a filament was 
found wrapping itself around the luminous dwarf VCC 165, 
which itself is probably orbiting around NGC 4216. The 
alignment of three dwarfs along an apparently single stream 
is also intriguing. Projection effects might be an explanation 
unless some of them are second-generation objects formed in- 
situ in tidal debris, or if all the dwarfs were accreted from the 
same cosmological filament. 

The dominant spiral itself does not show any strong 
evidence of having experienced a recent merger or being 
currently involved in a tidal interaction with a massive 
companion. The stellar streams most likely result from 
the disruption of dwarf satellites, some of which are still 
visible but look highly highly disrupted on the NGVS images. 

While dwarf galaxy accretion, with an associated high rate 



of satellite destruction and the formation of stellar streams, is 
a priori expected for massive galaxies - cosmological semi- 
analytical models and simulations predi ct that they grow 
at low redshift through minor mergers (iNaab et al.l l2009l 
iKhochfar et aP 1201 U iMcLure et al.1 12011) - their presence 
in a cluster environment is rather surprising. A possible 
explanation is that NGC 4216 is the central galaxy of a group 
that has recently fallen into the Virgo Cluster, and that the 
merging activity seen in its surroundings might be considered 
as evidence for pre-processing. This scenario does not fully 
explain the apparently truncated Hi distribution of the spiral 
galaxy, however, if caused by a ram-pressure stripping event. 
The level of ram-pressure required to account for the Hi 
deficiency might imply that the spiral has already passed 
closer by the central regions of the cluster at least ~^1 Gyr 
ago. This hypothesis would raise the question of the survival 
of stellar streams in the cluster environment, which might be 
longer than once believed. However, other processes that do 
not request an interaction with the central ICM are also able 
to truncate the Hi disk. Follow-up observations, in particular 
a spectroscopic survey to get dynamical information on the 
system, would be needed to determine how long ago the 
group of NGC 4216, including its numerous satellites, joined 
the Virgo Cluster 

We have carried out a preliminary comparison of the fine 
structure index of NGC 4216 with that determined in other 
spiral galaxies in the Virgo Cluster, but also around massive 
galaxies in the field and in groups. NGC 4216 seems to have a 
particularly high number of streams typical of minor mergers, 
and associated dwarfs caught in the act of being disrupted. 
We intend to confirm this as part of other on-going deep sur- 
veys with the Megacam camera on the CFHT. The signifi- 
cance and implications of our results should also be probed 
using numerical simulations of galaxy evolution which may 
tell whether NGC 4216 has indeed had an unusually rich re- 
cent accretion history, if on the contrary its fellow-members 
have had a much less important minor merger activity than 
predicted, or whether our knowledge of stellar streams and its 
associated time scales should be revisited. 
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